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9£ : ABSTRACT 

Os . 

A model is constructed for the mass and dynamics of M87 and the Virgo Cluster. 
q | Existing surface photometry of the galaxy, mass estimates from X-ray observations of 

^ | the hot intracluster gas, and the velocity dispersions of early-type Virgo galaxies, all are 

used to constrain the run of dark matter density over < r <^ 2 Mpc in the cluster. The 

"universal" halo advocated by Navarro, Frenk, & White, pdm oc (r/r s ) _1 (l + r/r s )~ 2 , 
■ provides an excellent description of the combined data, as does a Hernquist profile 

with pd m oc (r/r s ) _1 (l + r/r s )~ 3 . These models are favored over isothermal spheres, 
fSj , and their central structure is preferred to density cusps either much stronger or much 

weaker than r -1 . The galaxies and gas in the cluster trace its total mass distribution 
qq | (pgai oc p gas oc Pdm)j the galaxies' velocity ellipsoid is close to isotropic, and the gas 

temperature follows the virial temperature profile of the dark halo. The virial radius 
i-^ ■ and mass and the intracluster gas fraction of Virgo are evaluated. 

6 ■ 

Subject headings: dark matter — galaxies: clusters: general — galaxies: clusters: 
individual (Virgo) — galaxies: individual (M87) — intergalactic medium 

> ' 

H , 1. Introduction 

Some recent simulations of structure formation in CDM universes have suggested that dark 
matter halos generally form with a "universal" density profile, in which p^m oc (r/r s )~ l (l + r/r s )~ 2 
(Navarro, Frenk, &; White 1997, hereafter NFW). However, this claim has been questioned by 
other numerical work (e.g., Anninos & Norman 1996; Kravtsov et al. 1998); in particular, Moore 
et al. (1998) find instead that p<i m — > r~ lA as r — > in a set of highly resolved N-body halos. 
The little observational evidence that has been brought to bear on this issue is also ambiguous; 
e.g., the basic form of the NFW halo appears to apply in CNOC clusters of galaxies (Carlberg 
et al. 1997), but not in the dwarf galaxy DD0154 (Burkert & Silk 1997). This Letter tests the 
viability of the proposed universal halo, among other models for the dark matter distribution, in 
the nearby Virgo Cluster {D = 15 Mpc; Pierce et al. 1994). 
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The X-ray emission from Virgo, which is centered on the cD galaxy M87, has been used many 
times to constrain mass models for the cluster core (Nulsen & Bohringer 1995, and references 
therein). Here, an attempt is made to supplement this standard analysis with larger-scale 
constraints on Pdm(f) from the spatial distribution and dynamics of the galaxies in Virgo. To do 
so, the infalling late-type galaxies that are responsible for much of the cluster's irregular optical 
structure are excluded from analysis. Instead, reference will be made only to the early Hubble 
types in Virgo, which are spatially concentrated (about a point ~ 1° NW of M87) with a roughly 
Gaussian velocity distribution (centered near the velocity of M87), and which therefore suggest the 
existence of a relatively smooth and relaxed underlying mass distribution (see Binggeli, Tammann, 
k. Sandage 1987). In what follows, this is essentially adopted as a postulate. It is further 
assumed that Virgo is spherically symmetric and centered on M87. This is a reasonable first-order 
characterization of the X-ray structure around the galaxy, and will serve as a coarse-average 
approximation to the overall cluster structure. As part of this, the apparent offset between M87 
and the centroid of the early-type galaxy isopleths is taken not to reflect a large perturbation 
in the fundamental mass distribution, and is ignored. While imperfect, these simplifications 
ultimately admit a model for Virgo that is of interest both for its basic structure — the central 
scaling p<i m cx r _1 of NFW is indicated, and is favored clearly over the steeper cusps of either 
Moore et al. (1998) or a singular isothermal sphere — and for the relations it suggests between the 
galaxies, gas, and dark matter in this cluster. 



2. Data and Halo Models 

A significant contribution to the total mass on the smallest spatial scales in Virgo comes 
from the stars in M87 itself. With the stellar (core) mass-to-light ratio of this galaxy known (van 
der Marel 1991), and taken to be independent of radius, the mass density profile pM87(r) follows 
from the -B-band surface photometry of de Vaucouleurs &: Nieto (1978), which extends to R ~ 100 
kpc. (Here and throughout, three-dimensional radii are denoted by r, and projected radii by 
R.) Specifically, fitting ps with density models from the family discussed by Dehnen (1993) and 
Tremaine et al. (1994) yields the profile defined in Table |l] (for D = 15 Mpc, and with As = 0.09 
mag of extinction [Burstein & Heiles 1984] taken into account). Note that this model closely 
approximates a standard i? 1 / 4 law in projection (see Dehnen 1993). 

Moving outwards, density and temperature profiles for the hot gas within 10 ^ r <^ 200 
kpc of M87 have been obtained from ROSAT PSPC observations by Nulsen & Bohringer (1995, 
hereafter NB95). Their analysis assumes homogeneity, or a single-phase medium, which could be 
a potential source of error if a cooling flow has produced a multiphase gas in Virgo; however, 
this possibility appears to be of concern only at radii near the low end of the range in the 
ROSAT data (viz. r ^ 15 kpc: Canizares et al. 1982; Tsai 1994). Thus, use is made here of the 
model-independent gravitating mass profile, M tot (r), which NB95 derive by applying the standard 
assumption of hydrostatic equilibrium to their p ga s(^) and T{r) measurements. 
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Still further out, the (azimuthally averaged) surface density of dwarf ellipticals in Virgo is 
given by Binggeli et al. (1987), who find N gal (R) oc [l + (i?/410kpc) 2 ]- 1 for 100 kpc £ R £ 1.5 Mpc. 
(The distribution of E/SO systems, though somewhat more irregular in detail, is also consistent 
with this average profile.) The offset between the X-ray and optical centroids of Virgo is small 
enough — about half a core radius, and not much larger than M87 itself — that the same N ga \(R) 
applies to galaxy counts taken in rings centered on either point. Binggeli et al.'s parametrization 
then implies a three-dimensional n ga i(r) oc [1 + (r/410kpc) 2 ] _3/ ' 2 . The galaxy data do not 
require this constant-density core, however; they are equally well fit by projections of either 
a Hernquist (1990) profile, n ga \(r) oc (r/1.1 Mpc) _1 (l + r/1.1 Mpc) -3 , or an NFW profile, 
n ga,i oc (r/560kpc) -1 (l + r/560kpc)~ 2 (cf. Fig. | below). 

Finally, Girardi et al. (1996) use the velocities of 180 early-type Virgo galaxies to measure 
the aperture velocity dispersion, <7 ap — the r.m.s. line-of-sight velocity for all objects within a 
given distance from M87 — as a function of R from ~ 200 kpc to 2 Mpc. Essentially the same 
o" ap (i?) profile is obtained with circular apertures taken about the optical center of Virgo instead 
(cf. Fadda et al. 1996). 

Given these data, a parametric model for the dark matter halo may be specified as 
Pdm( r ) = K f(r/r s ), with / a dimensionless form function. It is then convenient to refer the 
halo's normalization and scale to the observed stellar quantities in M87 (Table |]). Thus, with 
D = ircpK/T bLb and r\ = r s /a, the mass interior to any radius is 

r/a \ 3 " 7 3 rl r - 



M tot (r) = T B L B 



1 + r/a 
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in which D and r\ are fixed by fitting to the observed mass profile at r ^ 200 kpc. If a good fit is 
achieved, no free parameters are left in the model, and it remains only to check its extrapolation 
to larger radii against the dynamics of the early-type galaxies. To this end, the observed n ga i(r) 
(from any of the functional fits mentioned above) and the model M to t(r) are used to solve the 
spherical Jeans equation for the radial velocity dispersion of the galaxies: 

2/3gal 



%al0-r( r ) = exp (- J 



dr 



r 



f°° GM tot ( f 2/3 gal \ 
J r n ga i— — exp (J — dx J dx 
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where /3 ga i(r) = 1 — o"|/o" 2 is the usual measure of orbital anisotropy. The dispersion along any 
single line of sight then follows from the projection 

N gal a 2 p (R) = 2 J r n gal a r 2 (r) (l - /3 gal -j J ^ _ R2 , (3) 

with N gal (R) = 2 J™ n gal (r 2 - R 2 )- l ' 2 r dr. Finally, 

rR \ / rR \ _1 



CT a P (^)= ( / N gaX al(x)xdx I \J N gaX {x)xdx) (4) 
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is compared with the observed aperture velocity dispersion profile at R ^ 200 kpc. If this 
comparison is unsatisfactory, but the basic assumption of dynamical equilibrium is retained, then 
the underlying form of /?dm(?") m ay be rejected. A potential complication is that the orbital 
anisotropy of the galaxies is not known a priori, and a given f{r/r s ) can safely be discarded only 
if it cannot reproduce the empirical <7 ap (ii) for any assumed (3 ga \(r). Fortunately, the Virgo data 
prove to be quite consistent with a simple f3 g3i \ = 0, i.e., isotropic orbits. 



3. Results 

The field of possible halo profiles can be narrowed somewhat by first requiring consistency 
with a small subset of the data: from NB95 and Girardi et al. (1996), 

Mtot(r £30 kpc) = (1.6t°;f) x 1O 12 M 
M tot (r£ 150 kpc) = (1.9 ±0.8) x 1O 13 M (5) 
o- ap (i? < 2Mpc) = 640±|j kms -1 , 

where all uncertainties are ~ 2-<r limits. Figure [l] illustrates tests of four different form 
functions against these three constraints only: in panel (a), f(r/r s ) = [1 + {r/r s ) 2 \~ l (an 
isothermal sphere with a core; cf. NB95); in (b), f(r/r s ) = (r/r s )~°' 5 (l + r/r s )~ 2 ' 5 ; in panel (c), 
f(r/r s ) = {r/r s )- 1 {l+r/r s )- 2 (the model ofNFW); and in (d), f(r/r s ) = {r /r s )~ 1 ^(l + r /vs)' 1 * 
(as in the N-body simulations of Moore et al. 1998). The different line types in each panel contain 
combinations of D and rj that are consistent (for /3 ga i = 0) with the 2-cr limits on the individual 
data points in equation @; the various dots identify models which yield the best estimates for 
both members of different data pairs. A minimally self-consistent model for p^ m must draw its 
normalization and scale length from the intersection of the three bands in the appropriate panel 
of Fig. |l], and ideally the three dots should coincide at a unique (77, D). 

Panel (d) of Fig. || shows that a central cusp as steep as r~ L5 (Moore et al. 1998) cannot 
easily account for the observed M to t inside both 30 kpc and 150 kpc: the implied M^ m oc r 
at small r is shallower than what is inferred from the observed excess of M tot (r) over the stellar 
Myi^{r). This is even more of a problem for the singular isothermal sphere (pdm oc r~ 2 and 
Md m oc r), which is altogether inconsistent with the X-ray masses (see also NB95). On the other 
hand, the mass constraints at r ^ 200 kpc are satisfied by any of the halos with central pdm oc r _1 
or shallower; but the weak cusps in panels (a) and (b) give the best agreement for scale radii 
r s that are so small [to limit a steeply rising Md m (r)] as to imply total halo masses which are 
too low to support the velocity dispersion of the Virgo galaxies on large spatial scales.0 Thus, 



2 This conclusion is insensitive to the assumed form of pdm as r — > 00. Also, while some shallow-cusp models with 
radial anisotropy (0 < /3 ga i < 1) can reproduce the observed <r ap at R = 2 Mpc, they are inconsistent with the full 
profile at smaller radii. 
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the "universal" halo of NFW, with its central pdm oc r and with /3 ga i = 0, emerges as the best 
candidate for a self-consistent description of the Virgo Cluster [panel (c)]. 

This simple analysis is limited in two respects. (1) At very small r <^ 10-15 kpc, where 
-^M87( r ) 3> Afjjmlr), the halo could depart from pd m oc r -1 and not conflict with any observations. 
And (2) little can be said about the behavior of pdm in the limit r — > oo; for example, the same 
procedure that leads to Fig. [l] shows that the profile of Hernquist (1990), pdm oc (r/r s )~ 1 (l+r/r s )~ 3 
(as suggested by the early simulations of Dubinski & Carlberg 1991), is just as acceptable as 
the NFW pdm ( r As) -1 (l + r/r s )~ 2 in Virgo. The best fits of these two models have different 
parameters (given in Table |l]), but show no real physical differences on ^ 2 Mpc scales. 

Figure || compares the full sets of data discussed in §2 with models using the best-fit NFW 
and Hernquist halos. The top panel shows the best X-ray estimate and 95% confidence limits for 
Mtot(r) (dashed line and bold curves), along with results from optical spectroscopy of the stars in 
M87 (Sargent et al. 1978; filled triangles) and from radial velocities of its globular clusters (Cohen 
& Ryzhov 1997; open triangles). These latter sets of data independently corroborate the models. 
(The globular cluster dynamics specifically will be discussed in a future paper.) The middle panel 
of Fig. |2] then compares the surface density profile of early-type galaxies with that of each fit to 
the dark matter halo. This shows — the analysis did not assume — that n ga i cx pdm in Virgo. This 
fact is used to compute aperture dispersion profiles that, for /3 ga i = 0, show excellent agreement 
with the data in the bottom of Fig. ||. It also allows the uncertainties on rj in Table |l] to be 
estimated by fitting the projection of Pdm( r ) to N gai \(R); the limits on D then follow from plots 
like panel (c) of Fig. ||. 

In either of these models, the virial radius of the dark matter halo — roughly, that within 
which its mean density is 200 times the critical density for closure — is r2oo = 1-55 ± 0.06 Mpc (for 
H = 70 km s" 1 Mpc -1 ), corresponding to a mass of M dm (r 20 o) = (4.2 ± 0.5) x 10 14 M . This 
emphasizes that the dark matter (and the gas) around M87 is associated — by construction — with 
the whole of Virgo; and this weakens previous arguments (Binggeli et al. 1987) which suggested 
that the cluster could not be relaxed. 

Finally, Fig. || shows that the hot gas within 200 kpc of M87 directly traces the distribution 
of dark matter (and galaxies) in Virgo: p gas /Pdm ^ 0.035 (L>/15 Mpc) 15 . If this is assumed to 
hold throughout the entire cluster, then the similarity of its total stellar and gas masses (Arnaud 
et al. 1992) implies a global baryon fraction of ~ 7%, roughly in line with other low-temperature 
clusters (Arnaud Sz Evrard 1998). The bottom of Fig. ||| also shows a close agreement between the 
observed gas temperatures and the virial temperature, kT(r) = 0.6m/f(j^(r), of the dark matter 
halo. Thus, simple virialization of gas that traces the dark matter in an r~ 1 halo suffices to 
explain the density and temperature structure of the intracluster medium at the center of Virgo. 
This suggests — if an NFW-type halo truly is universal, and if the proportionality p gas oc pd m is 
typical — that central density cusps and temperature drops in the gas of X-ray clusters may not, 
by themselves, prove the existence of cooling flows. That said, however, radiative cooling times are 
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unequivocally short in the core of Virgo (t C ooi & 10 10 yr for r ^ 70 kpc); and it remains to be seen 
whether a single-phase gas virialized in an NFW halo can account for the low-energy X-ray line 
fluxes observed there, or if these still require mass drop-out and the development of a multiphase 
medium at small radii in a cooling flow (cf. Tsai 1994). 

Of course, Virgo by itself cannot establish the "universality" of a central pdm oc r _1 structure, 
even among only cluster-sized halos. It is worth noting, then, that clusters at z ~ 0.3 in the CNOC 
survey are also well fit by an NFW model for the dark matter (Carlberg et al. 1997). Moreover, 
the NFW concentration of Virgo, c = r2oo/ r s = 2.8 ± 0.7, is consistent with that of the CNOC 
clusters. This c seems low, however, by comparison with A-CDM simulations of halo evolution 
(NFW; Thomas et al. 1998). The physical origin of this discrepancy — and, for that matter, of the 
r _1 cusp itself — is unclear. 

Marisa Girardi provided the velocity dispersion data plotted in Figure ||, and an anonymous 
referee offered many helpful comments. This work was supported by NASA through grant number 
HF-1097.01-97A awarded by the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., for NASA under contract NAS5-26555. 
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Fig. 1. — Tests of four different models for the dark matter halo in Virgo. In each panel, solid 
lines contain all (rj,D) pairs which (by eq. O) satisfy the 2-a limits on M tot (30 kpc) ; dashed lines 
contain those which satisfy the limits on M tot (150kpc); and dotted lines contain those which (by 
eqs. 0-0) satisfy the limits on fj ap (2Mpc). Circles in each case mark the single (r],D) pair that 
yields the best estimates for both Mt o t(30kpc) and Mt o t(150kpc) [there is no such pair for the steep 
central cusp in (d)]; triangles yield the best M tot (30kpc) and the best <r ap (2 Mpc); and squares give 
the best M tot (150kpc) and cr ap (2Mpc). 

Fig. 2. — Comparison of best-fit mass models with the Virgo data (scaled for an assumed distance of 
15 Mpc), for both NFW and Hernquist (1990) dark matter distributions. Parameters for the NFW 
halo follow from panel (c) of Fig. [l]; those for the Hernquist halo come from a similar procedure. 
The dotted lines in the top panel trace the separate mass components Mms7(0 an d Mdm( r )- in 
the middle panel, the profile of E/SO galaxies has been scaled to that of the dE's according to the 
ratio (697/56) of their total populations in Binggeli et al. (1987). 

Fig. 3. — Comparison of dark matter and intracluster gas density and temperature profiles. Data 
and (2-a) errobars are those from rings 2 through 37 of NB95. Points in the innermost annulus 
are somewhat uncertain, due to possible contamination by emission from the center of M87, while 
densities in the outermost one or two annuli may suffer from an inaccurate correction for background 
emission (NB95). Curves in the bottom panel trace the virial temperature profile of the dark 
matter halo, in the total potential arising from /?M87 + Pdm- See NB95 for discussion of the T(r) 
measurements. 
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Table 1. Mass Model for the Virgo Cluster 



PM87 



(i + S) 



7—4 



M M87 = T B L B {0^ 



3-7 



7 = 1.33 a = 5.1 ±0.6 kpc 
L B = (5.55 ± 0.50) x 10 10 L &)B 
T B = 14.6 ± 0.2 MqL~] b 



-i 



M dm = 4vr^3[ln(l + i r )-( T I^ 7 ) 



P^ = K(^y 1 (l + i 
M dm = 2vrKr3( T ^-)^ 



if 



= (3.2±f;|) x lO- 4 M pc~ 3 
= ira 3 K/T B L B = (1.65+ L35 ^ 



r s — 560 + 2oo 



\n./ = ^.uutalo) X 10 4 

i;^o kpc v = rs/a = no;30 



(1.6 ±0.7) x 10- 4 M Q p C - 3 
ira 3 K/T B L B = (8.1 ± 3.6) x 10~ 5 

1 ^-0.20 A,r__ . __ __ /. mn-4( 



if = 

D = 

r s = 1.07^g;§g Mpc r? = r s /a = 2107^ 
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